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Abstract 

The velocity anomaly recently reported by the OPERA collaboration appears strikingly at odds 
with the theory of special relativity. I offer a reinterpretation which removes this conflict, to wit 
that neutrinos yield a truer measurement of Einstein's limiting speed, and that light and indeed 
all other matter are retarded by additional interactions with the dark universe. I discuss existing 
experimental constraints and show that such a notion, considered cosmologically, can be subsumed 
in the dark-energy equation of state in an expanding Friedman- Robertson- Walker (FRW) universe. 
Planned measurements of the temporal variation in rcdshift have the potential to distinguish the 
possibilities. 



The OPERA collaboration has recently observed, with a significance of G.Ocr if statistical 
and systematic uncertainties are combined in quadrature, that neutrinos traverse a known 
length faster than they would were the speed of light in vacuum assumed [Tj, signalling 
an apparent violation of the theory of special relativity. In the face of this extraordinary 
outcome [2] alternate possibilities must be considered. Let us set aside obvious and admit- 
tedly more probable possibilities from the onset: we shall assume that the significance of the 
experimental result is robust and thus not obviated by unknown systematic errors, and that 
the OPERA collaboration has measured the group, rather than the phase, velocity of the 
"muon neutrino" wave packet. The empirical reality of neutrino oscillations [3J complicates 
the simple picture of a propagating z/^; the propagating mass eigenstates are not those of 
fiavor. Nevertheless, we assume that such details are not important here, so that the appar- 
ent violation of special relativity is manifest. In this context, then, is it possible to reconcile 
the OPERA result with the theory of special relativity? An affirmative answer requires an 
adjustment in our way of thinking. Consider that terrestrial measurements of the fundamen- 
tal speed of light are made under conditions in which only known matter is clearly absent, 
but cosmology tells us we live in a dark-energy and dark-matter dominated universe [:4j. I 
propose that neutrinos interact more weakly with the dark universe than photons and all 
other known matter do, so that propagating neutrinos offer a better measure of Einstein's 
limiting velocity. In this picture the photon remains massless, so that classical electrody- 
namics is unaltered, but interactions with the dark universe retard its speed slightly. This is 
tantamount to an index of refraction which differs slightly from unity, so that c can be less 
than the neutrino speed. This possibility is not at odds with special relativity, for the 
photon speed is measured in a background of dark energy and dark matter, rather than in 
a vacuum devoid of such content. We know little about the nature of dark matter, and less 
about dark energy; nevertheless, severe constraints do exist on the nature of the interactions 
we posit with the dark sector. We consider them carefully in what follows before turning to 
cosmological tests of the picture we espouse. 

Present Constraints. — The empirical value of c is 2.99792458 ■ 10® m/s [5]. We employ it 
as a conversion factor to relate time to distance, and indeed use it to define the meter. Its 
precise numerical value, however, does not derive from any known fundamental principle, 
so that it in itself is not sacrosanct. Stringent experimental tests, both terrestrial and 
cosmological, speak to the nature of the speed of light and hence to constraints on the 
interactions with the dark universe we propose. Severe limits exist on the variation of the 
speed of light with frequency [6l [7j, as well as on variations in the speed of light with respect 
to the orientation of its velocity vector in space [H [9] . Moreover, its universality as a limiting 
value of speed is also well-established [UlTO]. Thus we suppose the needed interaction must 
be energy-independent, isotropic, and universal for all matter save neutrinos. Under these 
conditions the speed of light, determined at the present cosmological time, can remain the 
same in every inertial reference frame presently accessible to us, albeit c < Cy. 

We assert that the speed of light c and the neutrino limiting speed Cy is related via 
c = Cjy/r;,, where n is an index of refraction with n > 1. In order to make our discussion 
concrete, we must employ an explicit framework for n. For definiteness suppose that we have 
a medium of scatterers, each of mass M, with mass density p. Employing the conventions 
and analysis of Ref. [11] we have, for a photon of angular frequency a;, 

for \n — 1| ^ 1, and where A^f is the forward Compton amplitude in the scatterer rest 
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frame. Assuming the discrete symmetries parity, time-reversal, charge-conjugation, as well as 
Lorentz invariance, unitarity, and analyticity, we have, for energies below particle-production 
threshold, that Mf = J2j=o^^j^^^ ^^^ere > for j > 1 pHj. To explain the OPERA 
data in the face of empirical constraints on the speed of light, we set all to zero save 
for j = 1, so that pA2/4M^ is set by their result for {v — c)/c, namely, pA2/4M^ = (2.48 ± 
0.28(stat)±0.30(sys)) x 10~^ = (5o [!]• We have posited a background of unknown matter, but 
our result, namely that the photon sees an index of refraction which differs from unity, can 
be of broader origin. In particular, the index formula can be generalized to particles of zero 
mass through the introduction of a thermal bath fl3[ [H] , and the structure of our expansion 
in u remains unaltered. The stringency of the tests on nonobservation of anisotropics in the 
speed of light suggest that our "unknown matter" is something other than the dark matter 
invoked to explain the observed galactic rotation curves. 

The OPERA experiment is not the first to study arrival time differences of photons and 
neutrinos. A previous short baseline experiment was sensitive to deviations in \v — c\/c 
to \v — c\/c < 4 X 10~^ [15], and the MINOS collaboration has reported a measurement 
of \v — c\/c = 5.1 ± 2.9 X 10~^ [16]. The observation of neutrinos from the supernova 
SN1987A [T7] sets a much more stringent limit. The neutrinos were observed to arrive 
some 3 hours before the first detection of optical brightening to yield a conservative limit 
of |f — c|/c ~ 2 X 10~^ [IE]. The limit implicitly assumes that the initial neutrino and 
photon pulses were emitted simultaneously, though we do expect the thermal neutrino burst 
from the core collapse, with neutrino energies of (9(10 MeV), to be emitted prior to the 
emission of visible light [19]. In Ref. [19] this time difference is assessed at ~ 10 hours, 
for a red-supergiant progenitor, though taking the OPERA result at face value implies 
that the neutrino burst associated with SN1987A was emitted some 4 years after optical 
emission. Such a long time lag is puzzling, and perhaps even implausible, but it must be 
noted that the detailed mechanism of a core-collapse supernova has not been established [20] . 
Moreover, SN 1987A in itself had many unusual features — e.g., its progenitor star was a 
blue supergiant [21]. The observed luminosity was also roughly an order of magnitude 
smaller than a typical Type II supernova and may be the result of the denser makeup of 
the star [2T]. Although the detected neutrino energies and burst duration appear consistent 
with their emission in core collapse leading to a proto-neutron-star [22] , the pulsar expected 
in this picture has not yet been observed. Alternate mechanisms are possible [23H2S]' 
features such as rotation of the core-collapse remnant and its associated magnetic fields may 
also play a role ^25j. In a black-hole-accretion-disk scenario, or "collapsar" model |23], e.g., 
neutrinos are emitted from the edge of the accretion disk formed after the core collapse; 
they are also of MeV energy scale — and the burst duration can be of comparable duration, 
though this outcome depends on the parameters of the model, as does whether the core- 
collapse neutrinos are trapped within the star [26l |27] . Particular features of SN 1987A 
suggest it may have had a companion star [28] as well; perhaps the dynamics of a binary 
system help explain the needed time lag and burst duration, with the close association of the 
observed neutrino and optical bursts attributable to coincidence. Observations of neutrinos 
from gamma-ray bursts could well yield more discriminating limits [2^] , but these have not 
yet been observed. 

The need to confirm the OPERA result in an independent experiment is clear. Beyond 
this, the empirical determination of an energy dependence in the limiting speed of the neu- 
trino would speak to complications beyond the simple picture we propose here. Recently 
constraints on superluminal models from the mere detection of neutrinos in the OPERA 
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FIG. 1: The change in recession velocity v as a. function of redshift z in various, flat ACDM 
cosmologies, as well as in some alternate scenarios. The solid lines from bottom to top denote 
= 0, 0.5, 0.729 (WMAP 7-year "best fit"), and 0.9, respectively, with Ho = 70.3 km/s/Mpcg2]. 
The dot-dashed lines bracket the associated solid ones above and below when w = — 1 is replaced 
by w = — 1.2 and w = —0.8, respectively. The dashed curve results when w is replaced by the 
z-dependent function used in Ref. [l3] with tdq = — 1 and Wa = 0.28; the dotted curves result when 
the calculation of v is amended by an index of refraction as per the OPERA result. 

experiment have been discussed j30]432] . These do not operate in the picture espoused 
here because Lorentz symmetry is not broken at the level of particle interactions. Pair 
bremsstrahlung, as discussed in Ref. [30], can nevertheless occur but via an explicit inter- 
action with the medium, much as computed in lepton "trident" production in quantum 
electrodynamics [33]. The observation of earth-crossing neutrinos in the OPERA energy 
range and beyond [3l] show that the presence of such pair production effects do not con- 
strain our scenario. Confirmation of our particular scenario requires study of the speed of 
light in temporally different regimes. To realize this, we turn to cosmological studies. 

Time Variation in Redshift. — The redshift to an object in a universe with matter and dark 
energy will change with time, and the measurement of its rate of change gives direct access to 
cosmological parameters [33HH7]. The measurement is very challenging, and may well require 
decades of observing time 1.36, '371, tiut with the advent of extremely large telescopes and 
laser-comb-stabilized calibration of spectrographs [38], it becomes possible [39]. It is useful 
to recap the standard computation of z [ID] before turning to the inclusion of interactions. 
In an expanding FRW universe, the redshift z of an object at time ti observed at time to 
is related to the cosmological scale factor a via 1 + z = a{to)/a{ti). Differentiating with 
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respect to to and noting dti/dto = 1/(1 + z), we find [36] 

dz 

z^ — = Ho{l + z)-H{z), (2) 

where the spectroscopic velocity shift is given hj i) = cz/{l + z) and H{z) is given in terms 
of 

H{z) = Ho [Qm{1 + zf + ^r{1 + zY + fiA(l + zf^'^""^ + (1 - aot)(l + zf] , (3) 

where Qm, ^r, and fl\ represent the fraction of energy density in matter, radiation, and 
dark energy, respectively, relative to the critical density today. For the region of z in which 
the stars and galaxies we observe reside, is completely negligible, and, moreover, we shall 
assume a flat cosmology so that Qtot = 1- The dark-energy density is characterized by an 
equation of state w, where w = —1 corresponds to the cosmological constant. The value of 
w need not be a constant in z, and in nonstandard cosmologies the scaling of the Qm term 
can also be modified, note Ref. [37] for a discussion. Curves illustrating the evolution of 
i) with z are shown in Fig. [TJ The possibility of light-dark-sector interactions can modify 
the shapes of these curves in z. The usual comoving distance is defined in the absence of 
interactions |41J; to include them via an index of refraction we note that the infinitesimal 
comoving distance is modified from 6t/a to n6t/a as the slowing of the speed of light makes 
the lightcone travel time longer [11]. Consequently our starting point is modified to 

a(ti) n(to) ^ ^ 

where we evaluate n{z{t)), noting in this case that n{z) = 1 -|- (1 + z)^6q [jTT]. Differentiating 
with respect to to and noting that citi/cito = ^/i^+z) as before, as well as that dz{tQ)/dtQ = 0, 
yields 



1 _ fl I ^) dlnniz) ^ ^ \ dz 

^ ' dz ' 

to leading order in small quantities, where v = Cyzj{\ -|- z\ though the replacement of the 
overall factor of c with in v will always be insignificant as it appears in a product with 
i^o- Note |n — 1| ~ (9(0.1) for z ~ 10. We compare the index-of-refraction-modified result 
for i} with that from the fiat ACDM model, as well as with other scenarios, in Fig. [Tj The 
inclusion of n shifts v to greater values with z^ just as a value of A > itself does [11]. We 
show, too, how V with z changes \i w ^ —1. Experimental constraints on w exist largely 
for models with constant w, so that, e.g., the WMAP 6-parameter ACDM fit to the 7-year 
data, combined with data from baryon-acoustic oscillations (BAO) and the measured value 
of iJo; yields w = —1.10 ± 0.14 [12]. In comparison, a direct measurement of the BAO 
angular scale using a distribution of galaxies with z = 0.5 — 0.6 yields w = —1.03 ± 0.16 if 
the other parameters are fixed j43j. This data set also yields a constraint on w{z): writing 
w{z) = Wo + Wa{l — I / {I + z)) and using the WMAP 7- year "best-fit" parameters [12] yields 
Wa = 0.06 ±0.22 [33]. Since the empirical data allow w < —1, we show how v changes if w is 
altered to w = —1.2 or w = —0.8 in Fig. [l| as well as if we employ w{z) and the WMAP 7- 
year "best-fit" parameters with Wa = 0.28 [43j. The effect of the modification with n on i; is 
rather small, though it begins to be appreciable for z in excess of z ~ 0.3 — 0.4. Interestingly 
this is within the window of Lyman-a forest studies, for which peculiar motions are known 
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to be negligibly small [39| HI]. If measurements of (9(lcm/s) can be made [38], then a 
decade of observations can resolve differences of (9(0.1 cm/s) — such a lengthy campaign is 
being planned [39j . 

Conclusions. — The sobering significance of the OPERA result pQ, coupled with decades 
of successful tests of special relativity [2], prompts us to consider alternatives which can 
be consistent with both. The experimental result c < need not reflect a breaking of 
Poincare invariance but, rather, could speak to light-dark-sector interactions which yield an 
index of refraction which differs from unity. The interactions must be energy-independent 
and isotropic and universal to all matter save neutrinos, to be consistent with existing 
experimental results. Nevertheless, the suggestion can be tested through the measured 
time- variation in redshift p5M37] : such studies permit the direct assessment of cosmological 
parameters. We have determined the amendment to z which appears in the presence of 
an index of refraction. To realize a simple but definite form of n{z) we have asserted 
that the photon couples to unknown matter, yielding n{z) ~ (1 + z)^. Reality could be 
much richer, and the possibilities include not only interactions with dark energy but also 
modifications to gravity itself. The interactions of which we speak may also be specific 
to very low redshift |l5]. Alternatively, perhaps c,y > c can mediate additional radiative 
effects in the very early universe [l6j. Moreover, the notion that a cosmologically local 
speed of light is tied to a dark energy model in which w{z) + 1 ^ is possible irrespective of 
whether OPERA is correct. The causal velocity itself could change with cosmological epoch, 
providing a context for the study of the time-dependence of other fundamental constants, 
such as a = e^/ he jl?]. Perhaps the OPERA result opens a new window on the dark universe 
— and the observational studies retain their interest even if it does not. The breadth of 
the possibilities underscores the importance of the observational measurement of both dark 
energy and its equation of state in a range of cosmological epochs. 
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